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Abstract 



We derive loop equations for the one-link correlators of gauge and scalar fields in the 
Kazakov-Migdal model. These equations determine the solution of the model in the 
large limit and are similar to analogous equations for the Hermitean two-matrix 
model. We give an explicit solution of the equations for the case of a Gaussian, 
quadratic potential. We also show how similar calculations in a non-Gaussian case 
reduce to purely algebraic equations. 
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1 Introduction 



Solving QCD in the limit of a large number of colors is a classic problem Recently, 
Kazakov and Migdal have noted that it may indeed be possible to find an analytic 
solution of the large limit of a certain version of induced SU{N) lattice gauge theory. 
In their model, the Yang-Mills action is absent at lattice distance scales and is induced in 
the continuum limit by the interaction of the gauge field with a heavy Hermitean matrix- 
valued scalar field. The eigenvalues of the scalar field behave as the classical "master field" 
in the large N limit. The possibility of solving this model has attracted wide attention 
and the task of unraveling its physical content is a subject of ongoing research. Besides 
being a candidate for QCD, it is an interesting example of a matrix model in dimensions 
greater than one for which a solution in the large N limit may be attainable. 

The Kazakov-Migdal model has action 



flld<P.ll[dU.y]e'^'^'''K (1.1) 

X xy 



where 

f/] = - E A^trV(0x) + E Nii(l),U,y(l)yUly , (1.2) 

X xy 

(j)x and Uxy (= U^x) N x N Hermitean and unitary matrices, which live on lattice 
sites X and links xy between neighboring sites, respectively. Here [dU] denotes the Haar 
measure for integration over U{N). Using the Itzykson-Zuber integral 

,,3) 

where a^i'e the eigenvalues of and x and A(0) = YliKji.'P^ ~ 4^) is the Vandermonde 

determinant, ( |1.1|) can be presented as the effective theory for a scalar field 



Z = / l[d(f)le^^«^'^^ (1.4) 

x 

with the action 

^eff[0] = -NY.Vi<Pl) + Y.lnA'i<l>x) + Y.^nI[(j>xAy] (1-5) 



xy 



(The Vandermonde determinants result from converting the integration measure for ma- 
trices to one for eigenvalues.) 

When is large this integral can be evaluated in the saddle-point approximation. 
The saddle-point equation is 

f = (1.6) 
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This model was formally solved by Migdal , who derived a singular nonlinear integral 
equation for the spectral density of the eigenvalues. This equation is rather complicated 
and its explicit solution even for the simplest Gaussian case proved to be non-trivial p. 
One of the authors has proposed an approach to this problem which employs the 
method of loop equations. In this paper we shall generalize these ideas to derive a loop 
equation for arbitrary one-link correlation functions. Its solution provides the spectral 
density of eigenvalues of cf) as well as the full set of one-link correlators of the gauge and 
scalar fields. We shall obtain an explicit solution of these equations in the case of a 
quadratic potential and reduce the more complicated case of a general potential to the 
solution of algebraic equations. 

The importance of the correlators of a pair of the gauge fields was discussed in |^ 
where it was argued that, to a large extent, they determine the continuum limit of the 
theory. More recently, in ^ it was shown that they could also play an essential role in 
resolving the problem of Z^q gauge invariance of the lattice model ( |1.1| ) |T0[ . Consider 



the correlator (our normalization differs from the original normalization of Ref. []^ by an 
extra 1 /N) 

Here, the delta-functions on the right-hand side arise from the gauge invariance of the 
integral. It also follows from gauge invariance that Cij is symmetric. Unitarity of U 
implies the sum rule 

^E^^. = l ■ (1-8) 

i 

In the large limit, the integration over the scalar field is performed by replacing it 
with its saddle point value (this should properly be done after is diagonalized) . Then, 
the expectation value of a quantity such as UW factorize into one-link expectation values. 

From this equation we can derive the quantity 

^T.C^AJ = ^^l^^['^'^]lx='^' (1-10) 

which can be obtained from the saddle-point equation ( |1.6|) as 

^Ec^.-0. = ^|m)-^E7x^l • (1-11) 




Formal expressions for the integrals ( |1.7| ) at finite A^, in terms of the eigenvalues 



have been found in [O. Those expressions are potentially very useful when N is small 
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rather than large. Here, to analyze the infinite limit, we shall take a different approach 
which requires a simultaneous solution for the correlator and the eigenvalue distribution. 

In the large limit, it is necessary to replace the index of the eigenvalues of (j) with 
a continuous label. This is done by introducing the eigenvalue density, p{a) such that 
p{a.)da is the number of eigenvalues in the interval [a, a + da]. The spectral density 
typically has support in a finite interval, [a, 6], and is normalized so that 

b 

dap{a) = 1 , (1.12) 



Also, the normalization of the correlator (|1.8| ) becomes 

rb 



dap{a)C{a,(3) = 1. (1.13) 



The saddle-point equation ( |1.11 



IS 



jyPp{P)C{a,P)P = ^ iV'ia) - 27^E„) (1.14) 

where TZE^ is the real (and continuous across the support of p) part of the analytic 
function (of the complex variable A) 



Equations of motion ( |1 . 1 1|) , (|1 . 14|) imply that, anywhere that the single-link operator 
Uxy4>yUyx appears with other operators inside a gauge invariant correlation function, and 
Uxy and (py appear no-where else in the other operators in the correlators, in the saddle- 
point approximation we can replace the diagonal elements of the operator Uxy(l)yUyx by 
the operator 

where V'{x) is the analytic continuation of 2TZE\ from the cut on the whole complex 
plane, so that on the cut 

V'{a) = 2nEo, . (1.17) 

This looks like the saddle point equation for the eigenvalue distribution of a one-matrix 
model with potential V{(j)) (see also Eq. (|2.16|) ). 



2 Loop equations for one-link correlators 
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2.1 The two analytic functions 

The central quantities of interest in this Section are the two analytic functions: Ex, which 
is defined by Eq. ( |1.15| ), and 

It is important to note that, due to invariance of the group measure under transformations 
U —>■ U\ Gy\ is symmetric in v and A. 

These two functions are analytic on the whole complex plane excluding cut singularities 
in the real axis in interval [a, 6] (the support of p). At the cut they have following 
discontinuities ^ 

Ea±io = V I dPp{P)^— T iT^pia) , (2.2) 
Ja a — p 

and 

G^i.cA = V f'dPpiP)^^ T t7Tp{a)gx{a) , (2.3) 
Ja a — p 

where 

g,{a)^ [dppiP)^^ (2.4) 

Note that this defines a function of the matrix (p^-, which we can also view as a 1-link 
expectation value 

Gx{<P.) = {Uxy^^Uyx)iL (2.5) 

A — (py 

(where the bracket denotes integration over Uxy and (py is evaluated at the saddle point). 
The function ^a(0) was introduced by Migdal ||^ and was the key component of his 
approach (see Sect. |^ below). 

Since the support of p is finite, the functions E^, and Gu\ have the following asymptotic 
expansions 

and 

As we see from Eq. ( pl7|) , asymptotics of G^x completely fixes Ex- 

2.2 Loop equation for Gj^x 

Our main equation is derived from the identity 

^ J d(l)[dU]tT-^e^ (^^j^Uxyj^Uyx) = {%-j^jf-j^Uxyj^Uyx) 

~ {wiT^^^^vT^^y^} ~^ {NUxz<f)zUzxi:z^Uxyjz^Uyx^ ■ (2.8) 
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The first term on the right-hand-side can be factored into two terms. This factorization 
is valid in the large hmit because of the saddle point evaluation of the integral over 0. 
The result is a product of two quantities 

/ tr 1 tr 1 1 \ 

The third term on the right-hand-side contains link operators which connect x to all 
neighboring points. One of these connects x to y. It contributes 

(^^—U^y^^Uy,) = -E, + AG.A (2.10) 

There are also 2D — 1 links which connect to other sites. For these one can use the fact 
that the quantity Uxz4>zUzx inside the expectation value bracket can be replaced by F^cj)^) 
(cf. Eq. ( P..16D ). With this input we present Eq. (|2.8|) in the form 

= E.G.A - (^J^L{<l>)^;^Uj^U^^ + XG,x - (2.11) 

Where, for brevity of notation, we have defined 

L{uj) = V'iio) - {2D - l)F{uj) = V'iio) + F{uj) . (2.12) 

If we assume that -^^(ci;) has no singularities in the interval [a,b], this equation can be 
written in the integral form ^ 

/ - — : — ^^^Gi^x = E,yG,yx + ^G^x — E,y , (2.13) 
Jci 2m u — 00 

where the contour Ci encircles the cut of the function Gi^x (which coincides with the 
support of the spectral density). 



2.3 Loop equation for E^, 

Consider the limit A ^ oo in Eq. ( ^.131 ). The leading term gives 

du L(uj) 



G.(.) = [ ^^±^E. - El (2.14) 

JCi 2m u — uj 



From the definition fp.lD of G 



\ N V — (p / \N V — (p / JCi 2-m v — lo 



•^An equivalent set of equations for an Hermitean two-matrix model was derived by Staudacher p2[ 
and by Ambjorn and Kristjansen [T^ . The point is that for the two-matrix model, which is associated 
with D = 1/2 in the formulas, one gets from Eq. ( 2.12D that L(iS) = V'iio). 
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Plugging this equation back to Eq. ( |2.14| ) and taking into account the definition ( |2.12|) 
one gets the equation for E^, first obtained in |T^ 

/ ^YMe^^eI (2.16) 

Jci 2m u ~ 00 

which looks like the large-A^ loop equation for an Hermitean one-matrix model with 
the potential V. However, in contrast to the one-matrix model where V is usually a 
polynomial, we shall consider below the case of non-polynomial V which has singularities 
on the complex plane outside the support of p. 



3 Relation to the Migdal's approach 

The purpose of this section is to demonstrate how, within our approach, one can recover 
the previous approach to solution of this model due to Migdal For simplicity we 
restrict ourselves to the case where the density of eigenvalues, p, has support in a single 
interval [a,b]. It is straightforward to generalize to the situation where the support is in 
two or more intervals. 

Consider our equation (|2.13| ). By construction (cf. Eq. ( p.lD ), the function G^x has 
singularities only on the support of the spectral density p, where it has a cut. The contour 
Ci encircles this cut, and the parameters u and A lie outside of it. Recall as well that the 
function L{uj) is assumed to be continuous across the cut (cf. Eq. ( p.l2| )). 



Let us compress the contour Ci to make it coincide with the cut. Then only the 



discontinuity of the function G^x at this cut (cf. Eq. ( |2.4| )) will contribute to the integral 
on the l.h.s. of Eq. (|2.13| ). Then Eq. ( p.l3| ) can be rewritten in the form 

y 



f dyp{y)L{y)^^ = E,G,x + XG^x - . (3.1) 

J a P — V 



Consider the limiting case of ( p.l|) when v approaches the cut singularity on the real 
axis from above: 

Vj',dyp{y)L{y)^ - znp{u)L{u)gx{u) = (3.2) 
{TIE, - trrpiu) + A) (v dyp{y)^ - mp{u)gx{y)) - UE, + iixp{u) . 



u-y 

Writing a similar equation in the limit when v approaches the cut from below and sub- 
tracting from it (|3.3|) , we finally get 

y-p 

Let us now make use of the relation between the real part of the function Ey and the 
function F{v) (cf. Eqs. ( |LT6D , ( [LTTQ , ( p32|) ) 

L{v) - F{v) = V'{v) = 2nEy , (3.4) 



V t dyp{y)^^ = _i + (A + TIE, - Liu)) Gxi^) (3.3) 

Ja V — U 
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where F{i>) is the logarithmic derivative of the Itzykson-Zuber integral which enters the 
saddle-point equation (|1.6|). Plugging ( p.4| ) into ( |3.3| ), we find exactly the basic equation 
of the Riemann-Hilbert approach of M 



V I dypiv)^^ = _i + (A - - 'JZE,)gx{u) . (3.5) 

Ja y — 
4 The exact solutions 



4.1 The general potential 

Let us consider an arbitrary potential V{uj), which is associated with a general ^^(u;) (for 
simplicity we consider the case when L{uj) is analytic on the whole complex plane, except, 
may be, infinity) 

oo 

Lico) = L^^"" . (4.1) 

m=0 

The equation for G ( |2.13| ) in this case is written explicitly 

{X + E,-L{u))G,x = E,-Rx{u) (4.2) 

where Rxiy) is given by 

Rx{^) = - / I^^G.A = E E ^^(A)^.™— 1 (4.3) 



m=l n=0 



and the contour C2 encircles both the cut singularity of and the pole ok, uj = v. 
The terms on the right-hand-side result from taking the residue at infinity in the contour 
integral and the functions G'n(A) are defined by Eq. ( p. 71) . 

The formal solution to Eq. ( [4.2|) is 

Gx= ^--^a(^) (44) 

with Rxiy) given by ([4.3|).P| The functions G'„(A) can be expressed in terms of Ex and 
L{ijj) using the recurrence relation 

G'n+i(A) = / ^.^^Gn{uj) - ExGniX) , Go(A) = Ex (4.5) 
JCi 2m A ~ 00 

which can easily be obtained by expanding Eq. ( |2.13| ) in 1/A. For n = this equation 
recovers Eq. (|2.14| ). 



It is worth mentioning that (4.4) obviously satisfies Eq. (pT5|) as is discussed in the previous section 



since Rx{v) has no discontinuity in v at the cut. 
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It still remains to determine Ex- It can be determined from the equation f\ 

^^^^L{u;)G^, = \E,-l (4.6) 



which is the l/u term of the asymptotic expansion of Eq. ( 2.13|) in i>. Using the expan- 
sion ( p.7|) , one can rewrite this equation as 

L^G^{\) = XEx-l. (4.7) 

m>0 

If L{X) were a polynomial of the highest power J, Eq. ( [4.7D and Eq. (|4.5| ) yield a 
polynomial equation for Ex containing powers of Ex up to order J. Some explicit solutions 
of this equation are given in the next subsections. Once Ex is known, V can be obtained 
from Eq. ( p^^ ). 

One could consider also the next terms of the l/i/-expansion of Eq. (|2.13| ) which result 
in the equations 

J2 LmGm+n{\) = AG„(A) + ^ EfcG„_i_fc(A) - E„ , n>l. (4.8) 

m>0 k=0 

A question arises as to whether these equations impose new restrictions on Ex- The 
answer is "no" and Eq. ( [4.8| ) is automatically satisfied as a consequence of Eqs. ( |4.5| ) and 
( [4.7[ ). The proof can be given by induction. Let us calculate the left-hand-side of Eq. 
at some n = tiq assuming that it holds for all lower n < Uq- One gets 



(4.5) (0) 

J2m>0 LmGm+n{X) = Z]m>0 -^m {/ci |^ (t^) — -E'AG'm+n-1 ( A) | = 

A {/c. tl^G^-ii^) - ExG^^^iX)} + YTuZl i?.G„_i_,(A) - j^. ^i^(^)G„_i(^) 
11) 

= AG„,(A) + Efc=o^fcG'„_i_fc(A)-E„ (4.9) 

which completes the proof. 

As is discussed above, G^x should be symmetric in v and A. While the right-hand-side 
of Eq. ( [4. 41 ) with Rx{i^) given by Eq. ( |4.3| ) looks non-symmetric, it can be shown that the 
symmetry is restored as a consequence of Eq. ( [4.71 ). To prove this, let us multiply the 
numerator and denominator on the right-hand-side of Eq. ( [4.4|) by {u + Ex — L{X))- The 
denominator then becomes explicitly symmetric and the numerator can be transformed 

■''The equation which was solved in Ref. [Q by the ansatz F{uj) = uj^L^^{lj) is just the 1/A^-term of 



the expansion of Eq. (4.6) in 1/A and is, therefore, incomplete. This ansatz is not consistent, generally 



speaking, with the whole Eq. (4.6). 
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as 

(4.5) 

{E,-Rx{jy)){jy + Ex-L{X)) = E,E^ + uE, - E,L{\) 

(il 

E^E, + uE^ - E.L(A) - E.Liu) + E^=o LmGm{X) + /c, /c. I^SS^- 

= E^E. - E^m - E.Liu) + uE^ + \E,-l + /^,^ J^, ^^G.. (.4.10) 

This expression is manifestly symmetric in u and A, what is easy to see taking the imag- 
inary part, which completes the proof. 

As is shown below for the quadratic and quartic cases, the symmetry requirement can 
be used directly to determine Ex alternatively to Eq. ( |4.7] ). 



4.2 Quadratic potential 

Let us show how the algorithm of the previous subsection works for the simplest case of 
a Gaussian potential. In this case, we begin with the assumption that L^ = for m ^ 1 
so that L{X) is 

L(A) = LiA = A-^A, (4.11) 

where we have introduced the quantity A to establish the connection with the previous 
studied of the Gaussian model. 

In the Gaussian case the solution ( |4.4| ) takes the form 

(0) _ E^ - A-^Ex 

since -Ra(^) = A^^Ex- Requiring this expression to be symmetric w.r.t. A and u one fixes 
the function Ex completely. It turns out that Ex coincides with the function obtained 
from the semi-circle eigenvalue distribution of the Gaussian matrix model (which is also 
known to solve the Gaussian Kazakov-Migdal model |^), 



(0) /iA ///2A2 1 /i2p 

while fj, and A are related by 

l^ = A-^-A. (4.14) 

Therefore, the symmetry requirement fully determines the Gaussian solution. 

This allows to rewrite G^^l in the manifestly symmetric form 

^(0) _ A~^ExE, - A~\XEx + uE,) + jXE, + uEx) + /i . . 

x^-(^A + A^i)Xu + u^ + fi ■ ^^-^^^ 
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It is straightforward to check that this expression is not singular for all v and A. 

As for the equation ( |4.7| ), in the Gaussian case it involves only Gi(A) which is fixed 
by the recurrence relation ( [4 .51 ) to be 

Gf{X) = \{\E^-\)-El. (4.16) 
This gives the quadratic equation 



whose solution, which at large A behaves like 1/A, is given by Eqs. ([4.13| ), (|4.14|) . 
Since (Eq. ([27[5| )) 

G.M=/ |^^£„=/"<i.pW^ (4.18) 

where the last integral is along the cut, F[x) is easily found from the discontinuity of 
G\[x) across the cut. The result coincides with that of 

fW(x)=Ax= (4.19) 

/i + VyU^ + 4 

It is crucial that D does not enter explicitly in the relation between /i and A (cf. 
Eqs. ( |4.14] ), ( [4.15| )). The meaning is that the Kazakov-Migdal model with the quadratic 
potential reduces to the Hermitean one-matrix model with the same kind of potential. 

Taking discontinuities of ( [4.15| ) across the cuts according to the definition ( p.l| ) we 
obtain the following expression for the correlator of the gauge fields 

Its expansion as — >■ oo recovers the leading order of the large mass expansion |^ 

Q . = 1 + 0,0 . . (4.21) 

On the contrary if ^ 0, one gets from ( [4.2(j| ) the weak coupling result |^ 

Ci,- = iV% . (4.22) 

While the expression (|4.20|) involves the quadratic polynomial in the denominator, its 
zeros, say in the complex A-plane, 



A,M = ^(A + -).±|^.^--. (4.23) 
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always lies for /i > outside the cut. 

The situation when C{a,l3) were become negative somewhere at the support of p, 
would be associated with a (third order) phase transition in the Itzykson-Zuber integral. 
The only possibility to have such a phase transition with our Gaussian Eq. ( |4.20| ) is when 
/i = so that \± = u belongs to the cut. Because M 



m^ = D^/i2 + 4-(Z}-l)/i, (4.24) 

the only possibility to reach the point /x = in the strong coupling phase mg > 2D, where 
the Gaussian model is stable, is at D = 1. This phase transition at D = 1 is the standard 
one which provides the continuum limit. 



4.3 The quart ic case 

One of simplest non-Gaussian cases is associated with 

L(A) = LiA + LsA^ . (4.25) 

The meaning of this formula is that we consider Eq. ( ^.25| ) as an ansatz for L{X) and shall 
solve Eq. ( |2.13D for Ex (and, therefore, V^'(A) = 27iEx at the cut). The original potential 
V of the Kazakov-Migdal model can then be determined due to Eq. ( |2.12| ) to be 

V\X) = 2DL{X) - {2D - 1)V'{X) . (4.26) 

As was already discussed, D = 1/2 for the two-matrix model and V disappears from this 
formula. For this reason one might think of / L{X)dX as of the potential of the proper 
two-matrix model which is just quartic for the ansatz ( |4.25| ). 



For L{X) given by Eq. (|4.25|) only the term with m = 3 survives in (|4.3| ) so that 
Eq. ( ^ gives 

^ Eu - (Li + L,u^)Ex - UuG.jX) - LsG2{X) 

Gux = TTT^ 7T~{ {^■2i) 

X + E^- L[u) 

where the recurrence relation (^.5|) yields 



Gi(A) = (Li + L^X^){XEx -I) -El- L-^E^ 
G2(A) = (L(A) - Ex) G'i(A) - XL-^^L^E^ + L^E, - 1) (4.28) 



with E2 and £"4 defined by Eq. (|27 

To determine Ex we use Eq. ( ^TT] ) (or, which is equivalent, the requirement that G^x 
should be symmetric in v and A) which leads to the following quartic equation 

Ex (A - LiL(A)) + L^G^iX) {Ex - L{X)) + L^El + L:,X\Lr + L3E2) = const . (4.29) 
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The constant on the right-hand-side of this equation can be found through asymptotic 
expansion in A of its left-hand-side. 

Since the resulting equation ( |4.29| ) for E,^ is quartic, its general solution, though explic- 
itly known, is rather obscure and not informative. For our purposes it is more convenient 
to rewrite it lowering the powers of Ex down to linear using Eq. which we rewrite 

in the form 

El = V\X)Ex - Q{\) (4.30) 

where 

du V'iuj) 



Jcv2 27rt A — uj 

is a polynomial of the highest power J — 2 if is \^ is that of the highest power J. It is 
essential for what follows that nether Q{X) nor V'{X) are singular at the cut of Ex. 

The resulting equation can be conveniently rewritten in terms of E{X) = L{X) — V'{X). 
Taking the discontinuity of the resulting equation {i.e. just the factor in front of Ex) one 
gets 



LiF(A) + L, [f^{X) - F{X)Q{X) - [Q{X) - {L, + L,X') - L3E2) (f(A) + L{X)] 
+XL^{L^E2 + - 1)] = A (4.32) 

while the continuous part yields 

Lig(A) + L3 [f{X)Q{X) (f(A) + L(A)) 
-(g(A) - (Li + LsA^) - L^E^) (g(A) - ^^(A)) 
-LsXL{X){LiE2 + L^Ei - 1)1 = L3A(L(A) + XL2E2) + const. (4.33) 



The set of equations ( [4.32|) , and ([4.33|) must have a solution with F{X) and Q{X) 



nonsingular at the cut. In general, Q{X) should be expressed via F{X), say, from Eq. ( |4.32|) . 
The substitution into Eq. (|4.33|) then yield a high order algebraic equation for -F(A). 
Presumably, a reasonable way to study the equation for -F(A) is to assume some analytic 
structure and then determine algebraic coefficients from the equation. From a bunch of 
solutions one should choose the one which recovers the Gaussian solution as L3 — > 0. 

While the general solution of Eqs. ( [4. 321 ), ( [4.33| ) looks hopeless, the iterative solution 
in L3 is straightforward 

and 

V'iX) = (^1 - + ^) A + ^3 (1 + A^ (4.35) 

Next orders in L3 can easily be calculated which illustrates that the solution exists at 
least at small enough L3. 
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Given the solution to Eqs. ( |4.32| ), ( |4.33| ), Cij can be calculated taking the double 
discontinuity of ( [4.27| ) w.r.t. both v and A across the cut which gives 



C{x,y) 



(Li + Lgx^) + L^xF{y) + L3 


'F\y) - Q{y) + {L, + Lsy') + L,E2 


y^-y(L{x)+F{x)) 


+ L{x)F{x) + Q{x) 



(4.36) 



As is discussed, this expression will be symmetric in x and y for F and Q being the 
solutions of Eqs. (|4.32| ), ([4.33|) . This requirement of the symmetry of (|4.36|) might help 
to find a solution. Anyway, the ansatz ( [4.25| ) was choosen as the simplest one while for 
other functions L{X) the resulting algebraic equations might simplify. 



5 Discussion 

In this Paper we have developed a method for calculating the gauge field correlator 
C{a,P). Its knowledge allows one to calculate the extended loop averages 

where UiT^y) is the ordered product of f/'s along some path T^y from x to y. If T^y 
coincides with one link, GuxiXxy) coincides with G^a according to the definition ( p.l|) . An 
explicit formula for GuxiT^y) can be obtained substituting the gauge field correlator at 
each link / by C{ai, a;+i) and integrating over a^'s. 

In the case of a matrix chain which is associated with the Kazakov-Migdal model at 
D < 1, this would be a complete set of observables at = cxo. However, a new class 
of correlators which can be constructed from pure gauge fields appears at D > 1. The 
simplest one is the closed adjoint Wilson loop 



Wa = 




(5.2) 



which is well-defined in the 't Hooft limit IQ and can again be calculated integrating over 
U's at each link. It is expected 0, that this quantity might undergo a first-order phase 
transition which is associated with a restoration of the area law. For the Gaussian model 
this phase transition does not happen while more complicated potentials have not yet 
been studied. 

One more interesting problem, for which the approach of this paper might be poten- 
tially useful, concerns another kind of pure gauge field averages, the filled Wilson loops 
which were employed in connection with the Kazakov-Migdal model in and whose 
importance for constructing the D > 1-dimensional string theories from the Kazakov- 
Migdal model was emphasized in and, as well as their role in resolving the problem of 
the local symmetry, in [R]. 
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We hope that the technique developed in this paper could help to solve some of these 
problems. 
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